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The surface roughness correction to the near-field heat transfer between two rough bulk materials 
is discussed by using second-order perturbation theory. The results allow for estimating the impact 
of surface roughness to the heat transfer in recent experiments between two plates and between 
a micro-sphere and a plate (using the Derjaguin approximation). Furthermore, we show that the 
proximity approximation for describing rough surfaces is valid for distances much smaller than the 
correlation length of the surface roughness even if the heat transfer is dominated by the coupling of 
surface modes. 
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I. INTRODUCTION 

Recently, several experimental setups made a measure- 
ment of the radiative heat flux on nanoscalc feasible^"— . 
It could be verified by Hu et al. that the heat flux be- 
tween two glass plates exceeds the far-field limit set by 
Planck's black body radiation lawi for distances of some 
microns. Narayanaswamy et al. and Shen et al. used 
a different experimental setup measuring the heat flux 
between microspheres and plates of different materials 
which allow for detecting the radiative heat flux at much 
smaller distances and have reported heat transfer coeffi- 
cients three orders of magnitude larger than the black 
body radiation limit^i^ in accordance with theoretical 
predictions^. A similar setup was used by Rousseau 
et al. for measuring the heat flux between a glass micro- 
sphere and a glass sample in a distance regime ranging 
from 30 nm to 2.5/zmi. With that experiment the theo- 
retical predictions based on fluctational electrodynamics^ 
could be verified with high accuracy. 

In all these experiments the results have been com- 
pared to calculations which do not take into account the 
roughness of the materials used. Here, we want to tackle 
the question how the surface roughness affects the heat 
flux in the near-field. Therewith we provide the basis 
for comparison of the experimental data with theoreti- 
cal results including roughness effects. Furthermore, our 
results might be used to study the impact of roughness 
in thermophotovoltaic devices^"—. The first work con- 
sidering surface roughness effects for the near-field heat 
fiux was given by Persson et ali^ employing the so-called 
proximity approximation (PA) ^"^d^ to determine the heat 
fiux between two rough surfaces. In a recent worki^ con- 
sidering the effect of roughness on the heat flux between 
a nanoparticle and a rough surface, it has been shown 
that perturbation theory exactly reproduces the PA for 
distances smaller than the correlation length of the rough 
surface. 

In this work, we extend the perturbation theory to 
describe the near-field heat transfer between two semi- 
infinite media with rough surfaces. This geometry is 
more suitable for the estimation of surface roughness ef- 



fects in recent experimental setups^. Utilizing the Der- 
jaguin approximation, our results can also be used to 
calculate the roughness correction for the experiments of 
Narayanaswamy et al., Shen et al^^ and Rousseau et 
alA. In addition, we show that the PA can be used for 
vacuum gaps smaller than the correlation length of the 
surface roughness even when the heat flux is dominantly 
due to the coupling of surface modes. This is an impor- 
tant result since it shows that the PA can indeed be used 
to obtain a simple estimation of the roughness correction 
to the heat fiux. 

The paper is organized as follows: In Sec. II we intro- 
duce the mean Poynting vector in terms of transmission 
coefficients and present in Sec. Ill the resulting expres- 
sions for second-order perturbation theory considering a 
Gaussian surface roughness. Finally, in Sec. IV we study 
the impact of surface roughness on the heat fiux and give 
a detailed discussion of the numerical results. In partic- 
ular, we investigate the range of validity of the PA. 



II. RADIATIVE HEAT TRANSFER 

Let us consider a configuration as depicted in Fig. [T] 
We have two semi-infinite bodies in local thermal equilib- 
rium at temperatures Ti and T2 separated by a vacuum 
gap. In general both media have different material prop- 
erties which can be expressed by different permittivies ei 
and £2- When considering isotropic and local materials 
only, the permittivities are scalars and do not depend on 
the wave vector. Here, we will consider such materials, 
although in general non-local effects have to be taken into 
accouniiiiiiS. 

Furthermore, we assume that both bodies have a rough 
surface expressed by the surface profile functions 5*1 (x) 
and S'2(x 4- d), where d is the mean distance between 
the two materials and x = {x,y). Since in this work we 
focus on stochastic Gaussian surface profilesi^, the profile 
functions fullfiU the properties 
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Figure 1: Sketch of the situation considered here. 

The brackets () stand for the average over an ensemble 
of reahzations of the surfaee profiles S'i(x) for i = 1, 2; 5^ 
is the rms height of the surface profiles. The correlation 
functions Wi(|x — x'|) are given by a Gaussian 



W,(|x-x'|) =e" 



(3) 



introducing the transverse correlation length a^. In addi- 
tion we assume that both surface profiles are statistically 
independent, i.e., {S1S2) = 0. For the Fourier compo- 
nent Si{K) of the surface profile functions one obtains 
for i = 1,2 

{S.{k))^0 (4) 
= (2^)''5?^(k + K')g,{K) (5) 
with the surface roughness power spectra 



(6) 



In order to evaluate the radiative heat transfer be- 
tween the two bodies^, it suffices to determine the electric 
Green's dyadic G(r,r') with r and r' within the vacuum 
gap. Then the mean Poynting vector is given by^ 



did 

2^ 



e(c^,ri)-e(w,r2) 



with 



{S^) = 2ReTr 



9,Gt(r,r')a,.G(r,r') 



(7) 



(8) 



where the integration is carried out over a flat surface 
within the vacuum gap parallel to the x-y plane and 
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(9) 



with I3i — \/{k-QTi) for t = 1,2; fee is Boltzmann's con- 
stant. Finally, we reformulate Eq. ([5]) in terms of trans- 
mission coefficients. To this end, we first introduce the 
Fourier representation of the Green's dyadic as 



G(r,r') 



d^K 



d^K' 
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with X = (x, y) and the lateral wave vector k = (fc^., ky). 
Then we find for (5^) in Eq. © 

(^")- E /(§/|;f^''-'^-^(^"'('^''^'-))' 

(11) 

where 

r„v(K,/«';a;) = 2ReTr <G{K.,K")d^d^><G\K' ,k") 



. (27r)2 , 

d,G\K,K")d,,<G{K\K" 



(12) 



is the transmission coefficient. It describes the trans- 
mission of a thermally emitted wave with wave vector k, 
and polarization j from body 1 to body 2 which is scat- 
tered into wave vector k' and polarization f during the 
transmission process. By means of the transmission co- 
efficient the Poynting vector in Eq. ([7]) can be converted 
into a Landauer-like form as was shown in Ref.— . Af- 
ter ensemble averaging, the translational and rotational 
symmetries are restored so that the mean transmission 
coefficient reduces to 



(T„.(/^,K';a.)) = {2iTf5{K-K')5,,,.T,{K-u:). 
and hence 
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In the following section we will determine the pertur- 
bation expansion of the transmission coefficient Tj = 

. and {S^} = sL''h{sL^^) + {sL^^)+. . . 



A. Perturbation result 

By using Eqs. p^ - p^ we obtain the perturbation ex- 
pansion for the mean Poynting vector by expanding the 
Green's dyadic with respect to the surface profile func- 
tions Si and 5*2. A short discussion of the perturba- 
tion method^ used here, can be found in Ref<2^. Note, 
that this perturbation approach has also been used to 
study for example the roughness effects for the Casimir 
force^^i^. We remark that the perturbation theory can 
be applied as far as the rms of the surface profiles is the 
smallest lenght scale of the problemiS^. 

Then, by inserting the zeroth-order dyadic Green's 
function from^ into the expression for the mean trans- 
mission coefficient in Eq. (|12p one retrieves for (Suj) from 
Eq. (fH)) the well-known result^ for the case of two media 
with flat surfaces 



fD~7P 

4Im(rj)Im(r|)c~^^'^ 
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fcg, fco = w/c and r' ( for z = 1, 2 and 



where 7^ = 

j = s, p) are Fresnel's reflection coefficients of the two 
bodies for s or p polarization; Djj = 1 — rjr| exp{2i'y^d) 



.2 



is the Fabry-Perot Uke denominator where 7^ 
tJ"-* is the transmission coefficient between the two bod- 
ies with flat surfaces having values between zero and one. 
In the propagating regime (k < fco) the property rj^-* < 1 
means that the heat flux is always smaller than that be- 
tween two black bodies for which Tj ^'^ = 1. In the evanes- 
cent regime the maximal value of one can be achieved for 
such pairs (w, k) which fuUfiU the dispersion relation of 
the coupled surface modes of the two surfaces^. 

To get the first nonvanishing term for the surface 
roughness correction, it is necessary to expand the 
Green's functions up to second order, since (G'^^) = 



and therefore T^^-* = (S'i'^'') = 0. The resulting second- 
order correction can be written as 
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where we have splitted the result into the specular part 
which depends on the mean field or mean Green's func- 
tion only, and the so-called diffuse part which is per def- 

(2\ (2'] 

inition given by {Sc^ ) — S'lLipcc- For the diffuse part, we 
obtainiSi 



2gi(|K - K'\){koSi)^ 



4|7,p|7;.p 

x[Q^^Qf{k-kr+Q'sQi'l^^ikxkr 



^2 l7l|' 



2 „2' I'f^'ci — • '^'7171 1^ 



+ (1^2), 



(17) 



where <*yp for i = 1,2 are the amplitude transmission 
coefficients of the two interfaces for s- and p-polarizcd 
modes; ^1 = it^ei — (i = 1,2), k = {kx,ky)/K and 

(1 o 2) symbolizes the term obtained by interchanging 

1/2 

the index 1 and 2. The functions Q , are defined as 

s/ p 



)l/2 
's/p 



s/p 
lAs/ppP 



+ 2Im(7,)Im(r^^/p')e-27d 



(18) 



It is seen from Eq. (jl7p that the diffuse correction to 

(2) 

the transmission coefficient and therefore also to ^Jjjg 
is always negative. Hence, due to the diffuse scattering 



in the rough surface the transmission coefficient becomes 
smaller and therefore the heat transfer less efficient. 

For the specular contribution to the mean transmission 
coefficient we obtain22 for propagating modes {k < fco) 



T, 



(2) 



(1 



+ (1 o 2) 



Im 



9 * - 

r2 e 



7r 



(19) 



and for evanescent modes (k > fco) 

(2 
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The expression for Ai/2,ss Ai/2,pp can be found in 

Refi^ and hV^ = 1,/ip^^ = £1/2- By comparing the 
specular second-order transmission coefficents with the 
zeroth-order expressions, it becomes apparent that they 



are very similar to t'^^^ , but with the difference that one 



of the terms 1 



.1/2: 



1 /2 

or 2Im(rj ) is replaced by a much 
more complex term involving Ai/2,jj- By means of this 
term the roughness scattering within one of the two sur- 
faces 1 or 2 is taken into account yielding a roughness 
correction to the transmission which can be either posi- 
tive or negative. 

In general, within second-order perturbation theory 
the diffuse and the specular components describe the cor- 
rection to the heat flux due to the scattering within only 
one of the rough surfaces, since both expressions split up 
into the sum of a term proportional to 51 and The 
different scattering processes can be classified as follows: 
The zeroth-order scattering is the specular scattering of 
the two flat mean surfaces. The diffuse component de- 
scribes the correction to the heat flux for first-order scat- 
tering of an incoming wave with lateral wave vector k into 
a wave with k' . For such scattering processes, incom- 
ing s-polarizcd waves can be scattered into p-polarized 
waves and propagating waves into evanescent waves and 
vice versa. On the other hand, the specular component 
describes the correction to the heat flux due to a second- 
order scattering of an incoming wave with k into a wave 
with k" ~ K, since it describes the scattering of the mean 
field. These scattering processes can be divided into two 

different types: a direct second-order scattering originat- 

~ f 2) 

ing from terms proportional to {k, — k") and an indi- 
rect second-order scattering through intermediate states 
with wave- vectors k' originating from terms proportional 
to sf'\K — K')sf'\K' — k"). The latter one is a sequence 
of two first-order scattering processes as for the diffuse 
scattering, but with the constraint that the polarization 
and the lateral wave vector are the same before and after 
the scattering sequence. 
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III. NUMERICAL RESULTS AND DISCUSSION 

In the following we discuss the impact of surface rough- 
ness on the near-field heat transfer numerically. We use 
the material properties of SiC and Au for both media 
with parameters taken from RefsJ^i^i and set T2 = OK 
and Ti = 300 K. Additionally, we choose the surface 
roughness parameters 6i = 5 nm and Oi = 200 nm for 
i = 1,2, i.e., we consider a very shallow surface rough- 
ness with 5i/ai ~ 0.025. 



A. Distance dependence 

First, we turn to the distance dependence of the rough- 
ness correction to the near-field heat transfer. To this 
aim, we start with a plot of the heat flux between two 
bodies with flat surfaces in Fig. [2] (a) considering two 
SiC and two Au slabs. SiC has a surface phonon reso- 
nance, whereas Au has no surface plasmon resonance in 
the infrared region. One can sec that the curves for the 
SiC and the Au plates arc quite different. For SiC the 
heat transfer is relatively large in the propagating regime 
(d > Ath = h/3c = 7.68 fj.m for T = 300 K) due to the low 
reflectivity of this material and has values close to the 
black body limit Sbb = 459.3 Wm~^, whereas for Au 
which has a high reflectivity the heat flux is very small. 
In the evanescent regime (d < Ath) the heat flux is for 
SiC dominated by the p-polarized surface mode contri- 
bution giving values which can exceed the black body 
result by more than 3 orders of magnitude for the here 
chosen distance regime. On the other hand, for the Au 
plates the heat flux is dominated by the s-polarized con- 
tribution and the large amount of transfcred heat can 
be attributed to induced Foucault currentsii. For dis- 
tances smaller than 100 nm the curve for SiC becomes 
proportional to which is the well-known result in 

the quasi-static limi t ^^1^^ . 
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t igure 2: (Color online) Plot of (a) 5^°' over distance for two 
SiC plates and two Au plates normalized to the black body 
value Sbb with (b) the corresponding roughness correction 
A{S). 

Now, we have a look on the impact of surface rough- 
ness. In Fig. [5] (b) we plot the roughness correction to 
the heat transfer AS* = (S'(2))/S'(o). First of aU, for 
both materials the roughness correction is positive for 
all distances so that roughness increases the heat flux 



with respect to the case of flat surfaces. In the propagat- 
ing regime the relative correction is very small (smaller 
than 0.01 percent) for SiC and is limited due to the fact 
that even if we introduce roughness the overall amount 
of thermal radiation cannot be larger than Sbb ■ For Au 
the relative correction is one order of magnitude larger, 
but the absolute correction is still very small, since 
has values smaller than 0.01 ■ Sbb- In the opposite limit 
of small distances, i.e., in the evanescent regime the rel- 
ative correction due to roughness can become large. For 
example, for SiC the relative correction is about 8 per- 
cent for c? = 40 nm. Since for that distance the heat flux 
is as large as 91 times the black body value, this is an 
increase of about 7 times Sbb- For the Au plates we 
flnd a roughness correction of 5 percent for d = 40 nm. 
At this distance 5''°-' is about 223 times the black body 
value yielding an absolute correction of about 11 times 
•Sbb- 
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Figure 3: (Color online) Plot of the specular and diffuse part 
of the roughness correction |A(S')| over distance for SiC. 

Finally, in Fig. [3] we plot the specular and diffuse part 
of the roughness correction for SiC only. Keeping in 
mind that the diffuse part is purely negative in this case, 
whereas we find that the specular part is purely positive. 
It can be seen that for distances larger than 500 nm the 
diffuse part is negligibly small compared to the specu- 
lar one, whereas for distances smaller than 500 nm both 
contributions have to be accounted for. Obviously, the 
specular part tends to overestimate the roughness correc- 
tion for such distances and only when also considering the 
diffuse contribution one will get the correct result. 

To understand why the diffuse part becomes important 
for small distances we take a closer look at the purely 
evanescent components of the diffuse Foynting vector in 
Eqs. (fT6|) and (fT7|) . For the scattering of evanescent p- 
modes into evanescent p-modes this component can be 
written as 

SL% = ' I d'^^ / A^K'Qlg2{\,^-K'\)f^^{K,K';uj)Q^^ 
+ {1^2) 

(21) 

where we have introduced the function /pp comprising all 
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the factors which are not essential for the discussion of 
the integral and 

Ql = r^^iM4>'"''- (22) 

The quantity is proportional to the transmission 

coefficient Tp"'' for evanescent p-modes coming from 
slab 1 which are transmitted into slab 2 and con- 
versely. This means that the integrand is proportional 

to Tp°''(k)p2(|'« — '^'\)tI'^\k') describing the transmis- 
sion of a p-polarized evanescent wave with wave vector k 
from slab 1 to slab 2 being scattered into a p-polarized 
evanescent wave with wave vector k' which is transmitted 
back into slab 1. Therefore, this product is proportional 
to the transmission coefficient describing how much of 
the energy is scattered back into the first slab by dif- 
fuse scattering within the second slab (see Fig. H]). In 
the non-retarded or quasi-static limit the weighting fac- 
tor 7exp(— 27^) for becomes KCxp{—2Kd), i.e., the 
transmission is the strongest if k w 1/d. Hence, one can 
conclude that the integrand is large when k w k' ^ 1/d. 
It follows that for — k'| ^ 1/d the integrand is small. 
Since the scattering within the rough surface is approxi- 
mately limited to ai|/« — k'| < 1 by the surface roughness 
power spectrum [see Eq. ([6])], we can conclude that for 

(2) 

tti <^ d the diffuse component ^jjj is comparatively 
small and can get large in the opposite limit, i.e., for 
d <^ ai. In addition, due to the Fabry-Pcrot denomina- 
tors in Qp and the diffuse contribution is especially 
large if k and k' corresponds to surface modes. Such 
scattering processes with k k' are reminiscent of direc- 
tional scattering of surface modes^i^. 




Figure 4: Sketch of the diffuse scattering in terms of trans- 
mission coefficients. Here Qp oc Tp^\k) and Q^' oc Tp'^\k'). 



B. Transmission coefficient in the (ixj, K)-plane 

To get some insight into the physical mechanism un- 
derlying the roughness correction, we will first discuss the 
transmission coefficients rj'"'' and rj'^'' which are defined 
in Eqs. (dH]) , (HZ]) , (HI]) and (HO]) and then we will turn to 



the spectrum of the zeroth-order Foynting vector Suj 

(2) 

and the second-order term {S^j ) ■ 

First, for distances larger or comparable to the ther- 
mal wavelength Ath ^ 7.68 /im at 300 K the propagating 
modes with k < ko dominate the contribution to the 
mean Foynting vector. For the two slab configuration 
these modes are given by the gap modes for which Tj^'^ 
equals 1 as shown in Fig. [5] (a) for p-polarized modes 
choosing d = 5 fj,m. The correction to the transmission 
coefficient due to roughness scattering is in this case dom- 
inated by the specular part. Nonetheless, this roughness 
correction is rather small for the chosen roughness pa- 
rameters as illustrated in Fig. [S] (b). 
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Figure 5: (Color online) Plot of (a) as defined in Eq. ^ 
and (b) Tp^cc/Tp"' in the propagating regime using d = 5 nm 
for frequencies ranging from cj = uJi = 1.83 ■ 10^** s~^ to 7 • 

Apart from this, due to the roughness scattering the 
evanescent surface polariton contribution can couple to 
the propagating modesi^ so that even for distances d > 
Ath one finds a relatively large correction to the transmis- 
sion coefficient for frequencies close to the surface res- 
onance frequency, which is for SiC given by wsphP = 
1.787 ■ 10^"* s~^. For d = 5 /im this correction is in the or- 
der of some percent. Since, the propagating modes domi- 
nate the heat transfer for d — b fim the overall correction 
to the Foynting vector is small. 

Now, for distances much smaller than the thermal 
wavelength Ath the near-field heat transfer is solely de- 
termined by the p-polarized surface mode contribution. 
For such distances the diffuse contribution cannot be ne- 
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glected. Choosing a distance of c? = 500 nm wc show 
in Fig. E] (a) a plot of Tp^^ first. Here, one can ob- 
serve the sphtting of the two surface phonon polariton 
branches^i^. Furthermore, one can estimate from Fig.[S] 
(a) that the relevant contributions stem from lateral wave 
vectors k smaller than about 1.5/a. In Fig. [5] (b) we 
show the corresponding plot of the roughness correction 



(a) 



Tp%:,cc/Tp^' . We find a qualitative similar correction as 
found for the roughness scattering considering only one 
rough surfacei^ with the difference (apart from the split- 
ting of the surface mode branches) that there is a positive 
correction for lateral wave vectors smaller than about 1/a 
for frequencies around wsphP and again a negative correc- 
tion for lateral wave vectors larger than 1/a. As we will 

(2) 

see below, when integrating Tp over the lateral wave 

vector to get (Sti ) this positive correction will make the 
negative one weaker compared to the case of only one 
rough surface in Refji&. This results in a positive overall 
correction to the heat flux for all distances in contrast to 
the results found for only one rough surface^. 

In Fig. [5] (c) we plot the corresponding diffuse contri- 
bution T^%g/T^°^ in the same (w, K)-range. As expected 
the diffuse component gives a negative correction which 
is relatively large for frequencies near the surface mode 
resonance. This correction can be as large as 2.5 percent 
for d = 500 nm. 

Hence, wc can conclude that in the evanescent regime 
where the heat transfer is due to the coupling of sur- 
face phonon polaritons the roughness correction can be 
relatively large even for a shallow surface roughnesses. 
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C. Spectral roughness correction 

Now, we turn to the effect of roughness scattering 
to the Poynting vector in the spectral domain, i.e., we 
concentrate on (S'if'')/S'w°\ In Fig. [7] we have plotted 

{s!^^)/s!;^^ for the specular and diffuse contribution us- 
ing d = 5 //m in a frequency range between 1.4 • 10^"* s~^ 

and 2 • lO^"* s^^ . We have also plotted S'i"'' using a scahng 
such that the curve fits into the plot. First of all one can 
see that; for such a large distance, the diffuse contribu- 
tion is negative as implied by Eq. ([T7| and about three 
orders of magnitude smaller than the specular contribu- 
tion. On the other hand, the specular contribution gives 
a relatively large correction of some percent near the sur- 
face resonance wsPhP = 1.787s~^, only. This correction 
is the sum of the positive and negative corrections to the 
propagating gap mode and the evanescent and propa- 
gating polariton mode correction within the reststrahlen 

region. Since the positive contribution to {s!if^)/s!j^ is 
larger than the negative one, one can expect to get a 
purely positive roughness correction to the Poynting vec- 
tor when carrying out the frequency integral. Further- 
more, since the heat flux is not solely determined by the 
surface mode, but rather dominated by the gap modes 
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Figure 6: (Color online) Plot of (a) T^°\ (b) T^%cc/T^°\ and 
(c) T^%ff/Tl'''' using d = 500 nm. 



(for frequencies greater than ui ) for which the roughness 
correction is very small, it follows that the roughness 
correction to the Poynting vector is small for distances 
around or greater than the thermal wavelength. 

In the evanescent regime, i.e., for a distance oi d ~ 
500 nm wc show in Fig. [5] the correction due to the scat- 
tering of surface phonon polaritons. As in the case for one 
rough surfacei^ we find for the specular part a positive 
correction for frequencies below and above the surface 
phonon polariton frequency (wsPhP — 1.787 s""'^) and a 
negative one for frequencies very close to the wsPhP- In 
contrast to the case of one rough surface, here the pos- 
itive contributions are larger than the negative ones so 
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Figure 7; (Color online) Plot of {SL^^)/S^^'' in a small fre- 
quency range around the surface phonon frequency of SiC 
(wsPhP = 1.787 • 10^* s~^) using d = 5/xm. We have plot- 
ted the specular (red solid line) and the diffuse (green dashed 
line) contributions (note the different scales at the right and 
left hand side) and s!^^ (blue dotted line). 
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D. Proximity approximation (d ^ min{ai, 02}) 

Now, we want to compare the numerical results for 
the heat flux with the so-called proximity approximation 
(PA). This will be useful to understand why the coherent 
flux is increased by the roughness. Within this approx- 
imation the rough surfaces are locally approximated by 
a flat surface (see Fig. ^ so that the mean flux can be 
computed from the flux between two flat surfaces S^^^ by 



(5(d))«(5W(d-5i(x) + 52(x))) 

= SpA- 



0(4) (23) 



When considering dielectric materials, then S^^^ oc 

in the quasi-static regime so that the PA can be further 

simplified to 



51) 



(24) 



that the overall Poynting vector is positive at this dis- 
tance. Furthermore, the diffuse contribution to the heat 
flux starts to play an important role and gives contribu- 
tions in the same order of magnitude as the specular one 
for frequencies near wsphp. This negative correction due 
to the diffuse component is still so small that the integral 
over all frequencies remains positive. Therefore, for the 
configuration considered here, the scattering of surface 
phonon polaritons docs not imply a negative correction 
to the heat flux as it was found for the heat transfer be- 
tween a nanoparticle and a rough surfacei^. In addition, 
since the zeroth-order sl^^ is dominated by the surface 
mode contribution for distances much smaller than Ath, 
the roughness correction of several percent gives a change 
of the heat flux of several percent. Hence, in the evanes- 
cent regime one can expect that the roughness correc- 
tion will be relatively large compared to the propagating 
regime. 
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For the case of two rough surfaces the PA was first used 
to estimate the near-field radiative heat transfer between 
two rough surfaces in Refii^. 
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Figure 8: (Color online) As Fig. [7] using d = 500 nm. 



Figure 9: (Color online) Sketch of the proximity approxima- 
tion. 

Now, with the perturbation result we can discuss the 
range of validity of this very simple approximation. As 
for the case of one rough surface discussed in Ref one 
can expect that the PA is valid as far as d <^ minjai, 02}, 
since then the local approximation by flat surfaces with 
a lateral extension much smaller than Oi can be assumed 
to be useful. In order to explore the range of validity 
of the PA we show in Fig. [TU]a plot of S{d)/SpAid) over 
the correlation length ai = a2 ~ a for different distances. 
It can be seen that the ratio S{d)/SpA{d) goes approxi- 
mately to 1 for a > d. Hence, we can conclude that the 
PA can be used for d < a. 

This conclusion is not trivial, since for SiC the near- 
fleld heat transfer is due to the coupling of thermally 
excited surface modes. The thermal near fields associ- 
ated with these surface modes can have a lateral coher- 
ence length Zcoh larger than the vacuum wave length^iS 
and can therefore be larger than the correlation length 
a considered here. However, in the quasi-static regime 
{d ^ Ath) it is found^i^ that Zcoh ~ d. Hence, for 
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Figure 10: (Color online) Plot of S{d)/S-PA{d) for SiC 
over a keeping the rms constant aA, 5 = 5nm for d = 
200 nm, 250 nm, 500 nm, and 750 nm. 



SiC converges for very small distances to the quasi-static 
expression for the PA in Eq. (j24p . since S'^^^ oc l/d? in 
the quasi-static regime. For Au 5^°^ saturates for small 
distances so that the quasi-static expression in Eq. ([M)) 
is not valid in that case. 

To summarize, for distances d <C min{ai,a2} the PA 
gives a good approximation of the impact of surface 
roughness for all the materials considered here. This al- 
lows for estimating the effect of surface roughness from 
the results obtained for a flat surface and is therefore 
an important result of this work. Note, that the same 
range of validity for the PA has also been reported for the 
Casimir forc e^'^i'^^'^^ . In this case, it has also been demon- 
strated that the theoretical predictions using the PA are 
in good agreement with the experimental dat a^^i"^^ . 



IV. CONCLUSION 
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Figure 11: (Color online) As Fig. [2] (b) but with the PA for 
SiC and Au calculated from with Eq. 

distances d <^ a the thermally generated fields above dif- 
ferent surface elements having an area smaller than a x a 
and larger than d x d are uncorrelated. Hence, for d <^ a 
the flat areas can indeed be regarded as independent so 
that the PA is valid. 

Numerically, we find for SiC that for distances d <^ a 
the correction to the heat flux is given by the sum of the 
positive specular and the negative diffuse contribution. 
The specular correction tends to give results larger than 
the one predicted by the PA. Only, when also considering 
the diffuse component of the heat flux, we retrieve the 
PA result. Therefore, the interplay of the diffuse and the 
specular scattering of the surface polaritons guarantees 
the validity of the PA. On the other hand, for Au which 
has no surface resonance around Ath we find that the 
diffuse part is negligible small. In this case the PA is also 
valid and is due to the specular part only. 

In Fig. [11] we plot the roughness correction for SiC and 
Au together with a numerical evaluation of the PA in 
Eq. ([23|. In both cases the PA gives a good approxima- 
tion for distances d a. As can be expected the curve for 



We have presented a perturbation theory for deter- 
mining the surface roughness correction to the heat fiux 
between two semi-infinite bodies. In particular, we have 
given the explicit expression for the mean Poynting vec- 
tor up to second-order perturbation theory assuming a 
Gaussian surface roughness. These results can be em- 
ployed to estimate the infiuence of roughness in recent 
experiments^"— and to study the impact of roughness in 
thermophotovoltaic devices^r— . 

In addition, we have discussed the numerical results for 
two SiC slabs with a given surface roughness. We found 
that the heat flux becomes larger when surface roughness 
is taken into account. By giving a detailed discussion 
of the correction to the transmission coefficients and the 
spectral Poynting vector for two SiC slabs, we could show 
that the scattering of surface modes within the rough 
surfaces causes a larger heat flux which is in contrast to 
the case of only one rough surface as discussed in Ref J^. 

Finally, we have shown that the PA is valid for dis- 
tances smaller than the correlation length of the surface 
roughness regardless of the materials (SiC, Au) consid- 
ered in this work. In particular, for the case of two SiC 
slabs the validity of the PA is due to an interplay of the 
specular and diffuse contribution to the mean Poynting 
vector, whereas for Au which do not have a surface po- 
lariton resonance in the infrared region the specular part 
already gives the PA result. Hence, the PA is valid even 
when the heat fiux is mainly due to the surface polariton 
coupling giving a simple and therefore powerful way of 
estimating the impact of surface roughness to the heat 
fiux. 
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